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ABSTRACT 

Using semi-analytical, one-dimensional models, we elucidate the influence of scattering and absorption on 
the degree of Ohmic dissipation in hot Jovian atmospheres. With the assumption of Saha equilibrium, the 
variation in temperature is the main driver of the variations in the electrical conductivity, induced current and 
Ohmic power dissipated. Atmospheres possessing temperature inversions tend to dissipate most of the Ohmic 
power superficially, at high altitudes, whereas those without temperature inversions are capable of greater 
dissipation deeper down. Scattering in the optical range of wavelengths tends to cool the lower atmosphere, 
thus reducing the degree of dissipation at depth. Purely absorbing cloud decks (in the infrared), of a finite extent 
in height, allow for localized reductions in dissipation and may reverse a temperature inversion if they are dense 
and thick enough, thus greatly enhancing the dissipation at depth. If Ohmic dissipation is the mechanism for 
inflating hot Jupiters, then variations in the atmospheric opacity (which may be interpreted as arising from 
variations in metallicity and cloud/haze properties) and magnetic field strength naturally produce a scatter in 
the measured radii at a given strength of irradiation. Future work will determine if these effects are dominant 
over evolutionary effects, which also contribute a scatter to the measured radii. 

Subject headings: planets and satellites: atmospheres 



1. INTRODUCTION 

A widely-accepted belief circulating within the astro- 
physical literature is that the intense starlight imping- 
ing upon hot Jovian atmospheres liberates electrons from 
their Group I metals (potassium and sodium), thus pro- 
viding a source of free electrons within the atmospheric 
flow (jBatygin & Stevenson! 1201 (A iPerna. Men ou & R auscherl 
120101) . The balance of collisional ionization with recombi- 
nation allows for the application of the Saha equation to- 
wards estimating the ionization fraction of the flow. If an 
ambient magnetic field exists, its intersection by the ~ 1 km 
s _1 flow produces induced currents which act to oppose the 
flow itself — a global, exoplanetary-scale manifestation of 
Lenz's law. The mostly neutral flow is coupled to the elec- 
trons via collisions, and the opposing forces act to retard the 
flow, a process known as magnetic drag. Ohmic dissipation 
converts the mechanical energy of the flow into heat, which 
acts to delay or even suspend the contraction of the exoplanet 
dWu & Lifhwickl2012l) . Thus, hot Jovian atmospheres behave 
like giant electrical circuits. 

It has been observed that hot Jupiters have radii which 
increase with the strength of irradiation, with a non- 
negligible scatter presen t in the measured rad i i at a given 
strength of irradiation ([Burrows et al.1 l2007ab lEnochetalJ 
201 1[ iLaiighlin. Crismani & Adamsll201 U iMiller & Fortnevl 
201 it iDemorv & SeagerJ 1201 ll) . (See also Figure 8 
of IPerna. Heng & Pontl 120121 ) A previously overlooked 
aspect in theoretical considerations of Ohmic dissipa- 
tion is that the variations of atmospheric scattering 
and absorption account naturally for this observed scat- 
ter, including the absence or pre sence of a tempera- 
ture inversion in the atmosphere dBurrows et al.1 [2007b; 
iFortney et al.1 12008). These variations in opacity may be 
attributed to variations in metallicity and/or the properties 
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of clouds/hazes (e.g.. iBurrows. Heng & Na mpaisarn 1201 lb 
iHelling. Jardine & Moklerll201 lb . the presence of which has 
been inferred from transi t observations obtained using the 
Hubble Space Tele scope dPont et alJl2008T: ISing et al.1 1201 lb 
iGibson et alEoTg) . 

In this Letter, we demonstrate these statements by using 
simpl e, semi-analytical models of t emperature-pressure pro- 
files dGuillodl20ia iHeng et al.ll2012l) to compute the Ohmic 
dissipation for a variety of model atmospheres. With the bene- 
fit of hindsight gleaned f rom three-dimensional sim ulations of 
atmospheric circulation (IPerna. Heng & Pond2012l) . we adopt 
simplifying assumptions for the atmospheric dynamics and 
focus on the variations in temperature, since the exponen- 
tial dependence of the ionization fraction on temperature is 
the overwhelming factor in determining the degree of Ohmic 
dissipation as well as how deeply it penetrates into the atmo- 
sphere. We describe our results in $2] discuss their implica- 
tions in [J3]and summarize our conclusions in Sj4] 

2. METHODOLOGY 

We assume a purely gaseous atmosphere dominated by 
molecular hydrogen and possessing an ambient magnetic field 
of unspecified geometry with a strength B. 

2.1. Time Scales 

Magnetic drag and Ohmic dissipation are effective pro- 
cesses only if the electrons and hydrogen molecules are 
tightly coupled, such that when the ambient magnetic field 
acts on the former, the effect is also felt by the lat- 
ter. For typical parameter values, the collisional time scale 
dDraine. Roberge & Dalgarnol 1983b . 
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is the shortest one in the system, where p denotes the mass 
density and T the temperature. In order for the magnetic 
effects to be communicated throughout the flow, we re- 
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FIG . 1 . — Top panel: temperature-pressure profiles used in this study, where 
the ratio of shortwave to longwave opacities (70) is varied. The irradiation 
temperature of all of these profiles has been set to T; rI . = f 850 K. Bottom 
panel: variation with pressure/height of the electrical conductivity (bottom 
left panel) and Ohmic power dissipated (bottom right panel) for models with 
70 = 0.1-10. For the bottom right panel, the horizontal axis at the top is the 
Ohmic power normalized by the incident stellar luminosity. 



(-B/1 G) _1 is the time scale associated with the gyration of 
an electron, with mass m e and charge e, about a magnetic 
field line. (The speed of light is denoted by c.) This condition 
is trivially satisfied, 



B < 3.4 kG 
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If the gyration time scale is shorter than the collisional time 
scale, then the opposing forces acting on the electrons are not 
communicated to the neutrals, but this requires implausible 
values for the magnetic field strength. 

Magnetic field lines anchored to the deep interior of a hot 
Jupiter may be perturbed by the atmospheric flow, across a 
pressure scale height H = kT/mg, with the perturbations 
being communicated at the Alfven speed va = B/2(np) 1 / 2 . 
The assumption of a steady-state, ambient magnetic field is 
reasonable if the Alfven time scale (£a ~ H /va) is less than 
the advective time scale associated with the zonal (east-west) 
atmospheric winds (i a dv ~ R/c s where R is the radius of the 
hot Jupiter and c s is the sound speed), which results in a lower 
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FIG. 2. — Models with the same value of 70 = 0.5 but with different values 
of the Bond albedo A. Left panel: temperature-pressure profiles. Right panel: 
Ohmic power dissipated for B = 1 G. For the right panel, the horizontal axis 
at the top is the Ohmic power normalized by the incident stellar luminosity. 



bound for the magnetic field strength, 
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We note that in all of the preceding estimates, the 
mass density is p ~ mg/kTKo ~ 10~ 6 g cm -3 
(T/1500 K) _1 (k /0.01 cm 2 g -1 ) -1 , where m = 2m H is the 
mean molecular mass, mn is the mass of the hydrogen atom, 
g = 10 m s~ 2 is the surface gravity and k is the Boltzmann 
constant. The quantity k is the broadband absorption opacity 
in either the optical or the infrared range of wavelengths. 

2.2. Model Description 

The incident flux impinging upon the substellar point is 
•Fq = "SB^lrr' wnere ^sb is the Stefan-Boltzmann constant 
and the irradiation temperature is 



Tin- = r* 
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Here, the stellar radius and effective temperature are rep- 
resented by R* and T*, respectively, while a is the spatial 
separation between the exoplanet and the star. We adopt 
Tin- = 1850 K, which is a value midway between that of HD 
189733b and HD 209458b, corresponding to Fq « 6.6 x 10 s 
erg cm' 2 s" 1 (1 - A) 1/4 ; on Earth, F « 1370 W m" 2 is 
the solar constant (if the Bond albedo A is set to zero). The 
incident stellar luminosity is £ = ^R 2 T^\ for R = 10 10 
cm, we get £ ~ 4.2 x 10 29 erg s _1 . 

For the radiative transfer, we adopt the dual-band ap- 
proximation, which assumes that incident starlight and ther- 
mal emissio n from the hot Jupiter occur at distinct wave- 
lengths (e.g lGuillotll2010HHeng, Frierson & Phillippsll201 li 
iHeng etalJl2012l) . Atmospheric absorption is described by 
two broadband opacities: n$ in the optical and «l in the in- 
frared. Their ratio is denoted by 70 = ks/^l- Atmospheric 
scattering in the optical is described by the Bond albedo A, 
albeit in the two-stream, collimated beam approximation. We 
ignore the effect of collision-induced absorption, although we 
note that such an effect increases the temperature, and hence 
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the Ohmic dissipation , at de pth and may be included via the 
models of Heng et al] (1201 2l) as an additional parameter. We 
set kl = 0.01 cm 2 g _1 such that the infrared photosphere 
resides at ~ 0.1 bar. The broadband, optical opacity then fol- 
lows from the assumed value of 70, i.e., ks = 7o«l. 

Dynamically, the zonal (east-west) wind speed v is the 
dominant component, such that the induced current density 
is primarily in the meridional (north-south) direction, 
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where a e is th e electrical conductivity. We ado pt the ex- 
pression used in lPerna. Menou & Rauscherl ((2010) for <r e and 
set the potassium abundance to be solar (ok = 10 -7 ). The 
Ohmic power dissipated in each layer of the atmosphere, with 
a thickness of one pressure scale height H, is 
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with T = 7/5 denoting the adiabatic gas index for an 
ideal diatomic gas. We have assumed that the zonal veloc- 
ity is equal to the sound speed, i.e., v — c s where c s = 
(TkT/mg) 1 / 2 . We truncate our models at P = 10 bar, 
since simulations of atmospheric circulation have shown that 
strong zonal winds, with speeds of v ~ 1 km s _1 , pene- 
trate down to this depth fo r the strength of stellar irradia- 
tion we have chosen (e.g., IHeng. Menou & Phillipps 1201 It 
iHeng. Frierson& Phillipps 201 1; Perna . Heng & Pontll2012l) . 



Since the zonal wind speeds approach or exceed the sound 
speed in these simulations, our assumption of v = c s is a 
reasonable one. 

2.3. Results 

The top panel of Figure Q] shows the temperature-pressure 
profiles for model atmospheres with various values of the op- 
tical to infrared opacity ratio, 70 = 0.1-10. Atmospheres 
with 70 > 1 pos sess temperature inversions, while those with 
7o < 1 do not dHubenv. Burrows & Sud arskv 20031: iHansenl 
2008; Guillot 2010). Even with this broad range in 70 values, 
the variations in temperature are modest, differing by less than 
a factor of 2. By contrast, the electrical conductivity spans 
several orders of magnitude, as shown in the bottom left panel 
of FigureQ] due to its exponential dependence on temperature 
via the Boltzmann factor. Atmospheres with temperature in- 
versions tend to concentrate Ohmic power at high altitudes 
(< 0.01 bar), while those without temperature inversions tend 
to dissipate more power deeper down (bottom right panel of 
Figure[T]i. For our model with 70 = 0.1, the local Ohmic dis- 
sipation at P = 10 bar approaches 0. 1 % of the incident stellar 
luminosity for B = 1 G; this value scales with B 2 . For the 
most highly irradiated hot Jupiters, the B 2 scaling is expected 
to saturate due to the effects of magnetic drag, which act to 
halt the zonal winds and thus reduce the Ohmic dissipation 
(lMenoull2012l) . 

We next examine the effects of scattering in the optical 
range of wavelengths. The left panel of Figure [2] shows the 
temperature-pressure profiles for models with 70 = 0.5 and 
three different values of the Bond albedo: A = 0,0.5 and 0.9. 
The right panel of Figure [2] shows the corresponding Ohmic 
power dissipated, as a function of pressure, for B = 1 G. 
Shortwave scattering produces two effects. The first and most 
obvious effect is to reduce the strength of the stellar irradia- 
tion via the diminution of T; rr . The second effect is to cool 
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FIG. 3. — Ohmic power dissipated at P = 10 bar for models with 70 = 
0.1-10 and A = 0,0.5 and 0.9. The assumed magnetic field strength is 
B = 1G. 



the lower atmosphere and thereby reduce the Ohmic dissipa- 
ti on at depth. 

Guill ot & Showman] (120021) have suggested that power dis- 
sipated at ~ 10 bar may affect the evolution of a hot Jupiter 
and provide an explanation for its infl ate d radius. (See 
alsolBurrows, Sudarskv & Hubbardl2003l and lWu & Lithwickl 
In Figure |5] we show the Ohmic power dissipated at 
P = 10 bar for our models with 70 = 0.1-10 and A = 0, 0.5 
and 0.9, again assuming B = 1 G. In addition to introduc- 
ing a spread in the Ohmic power dissipated at a given 70 
value, shortwave scattering also reduces the difference in Pj 
between low- and high-70 models for a given value of A. 

Finally, we consider purely-absorbing cloud decks (in the 
infrared) with a Gaussian shape/profile (Figure @|, as de- 
scr ibed in IHeng et al] (2012). We examine dense (k Co =0.1 
cm 2 g _1 ) versus tenuous (k Co = 0.01 cm 2 g _1 ) and thick 
(A c = 1) versus thin (A c = 10) cloud decks. In all of the 
models shown, the photon deposition depth (in the absence 
of scattering) is located at Pd ~ 0.63<?/ks = 31.5 mbar. 
We have intentionally centered t he cloud decks at P c = 31.5 
mbar, as it has been shown by Heng etail (1201 2h that this 
maximizes the effects we are about to discuss. Cloud decks 
absorbing purely in the infrared produce a greenhouse ef- 
fect, warming the lower atmosphere and cooling the upper 
atmosphere, with denser clouds producing a greater effect. If 
the cloud deck is thin enough, it imprints its shape onto the 
temperature-pressure profile, thus producing localized reduc- 
tions in the Ohmic power dissipated. If the cloud deck is dense 
and thick enough (A c = 1), then it reverses the temperature 
inversion altogether, thus greatly enhancing the Ohmic dissi- 
pation at depth. If the cloud deck sits above the photon de- 
position depth (P c < Pd), then these effects are diminished. 
If it sits below the photon deposition depth (P c > Pd), then 
it cools the upper atmosphere but has no effect on the lower 
atmosphere. 

We note that while the deposition of heat in the lower 
atmosphere pushes the radiative-convective boundary to 
larger pressures, thus preventing the exoplanet from los- 
ing i ts primordial entropy quickly dGuillot & Showman I 
2002|, this effect can only explain observed radii up to 
about 1.3 Jupiter radii. To explain the existence of 
hot Jupiters with larger radii, the dissipati on of heat in 
the deep, convective interior is required dBurrows et alj 
I2007at iBatvgin. Stevenson & Bodenheimerl 1201 lh . Thus, 
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FIG. 4. — Models with a Gaussian cloud deck, absorbing only in the in- 
frared, centered at P c = 0.01 bar. Top panel: temperature-pressure pro- 
files which include cloud decks with different absorption opacities and thick- 
nesses. Bottom panel: Ohmic power dissipated for B = 1 G. In both panels, 
the thin, yellow curve shows the cloud-free case for 70 = 2. 

even temperature-pressure profiles which dissipate most of 
the heat superficially (70 > 1) will still act to keep the ex- 
oplanet inflated because the deep isotherm corresponds to a 
non-negligible amount of electrical conductivity. 

We conclude that the variations in the absorption opacities 
(70) and albedo (.4), as well as the cloud/haze properties, nat- 
urally account for a range of values of the Ohmic power dis- 
sipated at depth (~ 10 bar). If Ohmic dissipation is the dom- 
inant mechanism for inflating hot Jupiters, then this implies 
a scatter in the measured radii at a given strength of stellar 
irradiation. 

3. DISCUSSION 

3.1. Connections to Previous Work 

We have developed models which elucidate the link be- 
tween atmospheric scattering, absorption and Ohmic dissipa- 
tion in hot Jupiters at a given strength of stellar irrad iation. 
Our study is complimentary to that of Menou (2012), who 
developed analytical scaling laws for the variation of mag- 
netic drag and Ohmic dissipation with the strength of stel- 
lar irradiation, while we have explored the effects of stellar 
irradiation on heat redistribution and energy dissipation (in- 
cluding Ohmic dissi pation), via numerical s imulations, in a 
companion paper (iPerna. Heng & Pontl l2012l). Our simplify- 
ing assumption for the atmospheric dynamics is corroborated 
by the results of these simulations. Furthermore, our work 



is also complimentary to that of IWu & Lithwickl d2012l) . who 
adopted a two-zone model consisting of an isothermal enve- 
lope (i.e., a constant temperature-pressure profile) surround- 
ing an isentropic core. The variety of temperature-pressure 
profiles used in the present study highlights the importance 
of atmospheric effects on the evolution of a hot Jupiter, as 
the atmosphere provide s a boundary condition for the evolu- 
tion of the interior (e.gjBurrows. S udarskv & Hubbar dl2003l ; 
Burrows, Heng & Nampaisarn 201 1). 

3.2. Evolutionary Effects 

We have not examined effects stemming from the evolu- 
tion of a hot Jupiter, which we will now discuss. Firstly, an 
increased atmospheric opacity slows down the gravitational 
contraction and internal cooling of the exo planet, as the flux 
eman a ted scales inversely wi th the opacity dArras & Bildstenl 
2006; Burrows et alJ l2007ah . Such a ret ardation or sus- 
pensio n of the cooling has been studied by IWu & Lithwickl 
d2012l) . who showed that Ohmic dissipation cannot effec- 
tively re-inflate a hot Jupiter after it has cooled. Secondly, 
an enhanced metallicity — whether in the form of a core 
or distributed throughout the envelope — implies an en- 
hanced surface gravity a nd a smaller radius of the exoplanet 
dMiller & Fortney 201 1). In turn, these effects imply that the 
amount of dissipation required to explain an observed radius 
may be orders of mag nitude larger compared to a situation 
with lower metallicit y ( Bodenheimer, Lin & Mardlin gjl2001l : 
iBodenheimer. Laughlin & Linl 120031 : iBurrows et alJ l2007al) . 
In other words, a hot Jupiter possessing a core or a metal- 
enriched atmosphere tends to be less inflated, all else being 
equal. Thirdly, the presence of a temperature inversion it- 
self will lead to smaller radii pu rely as an evolutionary ef- 
fect dParmentier & Guillotll201 lb . All of these evolutionary 
effects will contribute to the scatter in the measured radii of 
hot Jupiters at a given strength of stellar irradiation. 

3.3. Other Observational Consequences 

The variations in atmospheric opacity are manifested in 
other observable ways. At a basic level, hot Jovian atmo- 
spheres are characterized by two time scales, namely the 
radiative (t ra d) and advecti ye (t a dv ~ R/c s ) ti me scales 
(IShowman & Guillotl l200l ICowan & Agol 1201 lh. It is a 
gener ic statement that the radiative time scale dGoody & Yungl 
1989), 

increases with depth or pressure, such that the atmosphere is 
dominated by advection at large pressures. (The specific heat 
capacity at constant pressure is denoted by cp; its value is 
14550.4 J K- 1 kg- 1 for T = 7/5 and m = 2m H .) At the 
top of the atmosphere (< 1 mbar), we have i ra( j <C t a dv such 
that the atmosphere is predominantly radiative. If starlight 
is mostly deposited at the top of the atmosphere, then heat 
redistribution from the day- to the night-side hemisphere is 
inefficient; if it is deposited mostly at depth, then heat redis- 
tribution becomes efficient. 

In the presence of scattering, starlig ht is predominantl y ab- 
sorbed at the photon deposition depth (Heng et al. 2012b, 3 
0.63g (I- A 



1 + A 



(8) 



In the presence of scattering, the photon deposition depth always sits 
deeper than the optical photosphere. In the absence of scattering, the two 
layers are coincident. 
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FIG. 5. — Ratio of advective to radiative time scales, evaluated at the photon 
deposition depth, as a function of the optical to infrared opacity ratio 70 and 
for three values of the Bond albedo .4 = 0, 0.5 and 0.9. The broadband, 
infrared absorption opacity is fixed at kl = 0.01 cm 2 g — 1 . The atmosphere 
becomes more radiative (or less advective) as t a< j v /t ra< i increases. 

It is thus plausible to evaluate £ a dvArad at P = P\>. Figure 
[5] shows that £ a dvArad is in general an increasing function of 
70. For a fixed broadband, infrared absorption opacity kl, an 
increasing value of 70 corresponds to an increasing value of 
the broadband, optical absorption opacity ks and thus a de- 
creasing value of Pd- At a given value of 70, increasing the 
Bond albedo A also decreases Pq, since the overall penetra- 
tion depth of an optical stellar photon, in the vertical direction, 
decreases when scattering is present. 

As the photon deposition depth shifts to higher altitudes, 
it samples an increasingly radiative part of the atmosphere 
— advection becomes more sluggish and heat redistribution 
becomes less efficient. The angular shift of the peak, from 
the substellar point, of the photospheric infrared flux ema- 
nating from the hot Jupiter — known as the "hotspot off- 
set" — is a proxy for the efficiency of heat redistribution, 



with smaller shifts corresponding to less efficient redistribu- 
tion. It is thus natural to expect that — all other things be- 
ing equal (e.g., g, kl) — hot Jovian atmospheres possessing 
temperature inversions, at a given strength of stellar irradia- 
tion, will have less efficient heat redistribution and a higher 
day- to night-side flux contrast, an expectation which has 
be en corroborated by th e atmo spheric circulation simulations 
of iPerna. Heng & Pontl (|2012|) . While these simulations do 
not include scattering in the optical, we expect atmospheres 
with non-zero albedos to produce a similar effect. Ultimately, 
stellar irradiation remains the main driver of heat redistribu- 
tion and ene rgy dissipation with opacity effects playing a sec- 
ondary role (Pe rna. Heng & PontJl2012h . 



4. CONCLUSION 

Collectively, the opacity effects (absorption and scattering) 
we have discussed, which may be attributed to variations in 
the metallicity and/or cloud/haze properties, as well as ex- 
pected variations in the magnetic field strength, naturally in- 
troduce a scatter to the measured radii of hot Jupiters at a 
given strength of stellar irradiation Ohmic dissipation is the 
dominant mechanism for maintaining radius inflation. Fur- 
thermore, the evolutionary effects we discussed in £13.21 will 
also contribute a scatter to the measured radii and it remains 
to be determined which set of effects provides the dominant 
source of scatter. Future multi-wavelength measurements of 
thermal phase curves, from a larger sample of hot Jupiters, 
will disentangle these effects. 
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